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1. Introduction 

In this talk I will give an overview over recent results obtained within the Hamiltonian ap- 
proach to QCD in Coulomb gauge. I will focus on the extension of this approach to finite tem- 
peratures and to the inclusion of quarks. My main interest here will be the finite-temperature 
deconfinement phase transition and the spontaneous breaking of chiral symmetry. 

Going to Weyl gauge Aq = canonical quantization of Yang-Mills theory yields the following 
Hamiltonian 

# = ^|(n 2 +fi 2 ). (i.i) 

Here ft = 8/ iSA is the canonical momentum operator and B is the non-Abelian magnetic field. 
The Yang-Mills Hamiltonian (JOJ) is invariant under spatial (time-independent) gauge transforma- 
tions. A quantity of central interest in quantum field theory is the vacuum wave functional, by 
means of which all correlation functions can be evaluated. This quantity is obtained by solving the 
Schrodinger equation 

Hy = Ey (1.2) 

for the lowest energy eigenstate. There have been various attempts to solve the Yang-Mills Schro- 
dinger equation ( |L~2| ) directly for gauge invariant wave functionals, in particular in D = 2+ 1, 
see Ref. [|l|] and references therein. One can give arguments that the Yang-Mills vacuum wave 
functional (in D = 3 + 1) can be approximated in the low-energy regime by [0] 



y[A] = exp 



Fh 



(1.3) 



So far one has not succeeded in determining the Yang-Mills vacuum wave functional in a gauge 
invariant way. A much more convenient way is to fix the gauge and for the purpose of Hamiltonian 
Yang-Mills theory Coulomb gauge is, in particular, convenient. The prize one pays is that the 
Hamiltonian is more complicated. In Coulomb gauge the Yang-Mills Hamiltonian is given by [|3|] 

H = \J ' (j~ 1 fl ± Jti ± +S l )+Hc, (1.4) 

where / = Det(— Do) is the Faddeev-Popov determinant with D = a +gA, A = f ac0 A c being the 
covariant derivative in the adjoint representation. Furthermore, 

H C = l -Jj-mm = ^J j- { P {-bd)-\-d 2 )(-Dd)- i j P (i.5) 

is the Coulomb Hamiltonian, which arises from solving Gauss's law (which is a constraint to the 
wave functional to guarantee gauge invariance) for the longitudinal momentum operator n". In 



Eq. ( |1.5[ ), p = p g i + p m is the total color charge, which contains beside the charge of the Yang- 
Mills field p g i = —ATI also the charge of the matter fields p m . By resolving Gauss's law gauge 
invariance has been fully taken into account and in Coulomb gauge each functional of the transverse 
gluon field A L , dA L = is, in principle, a physical wave functional. Note that in the canonical 
quantization of Yang-Mills theory the gauge field figures as the coordinate, and the transition to 
Coulomb gauge implies a transition to curvilinear coordinates. Accordingly, the kinetic piece of 



2 



Hamilton Approach to QCD in Coulomb Gauge 



Hugo Reinhardt 



the Yang-Mills Hamiltonian [the first term in Eq. (1.4)] resembles the Laplacian in curvilinear 
coordinates. In the scalar product of the Yang-Mills wave functionals the transition to Coulomb 
gauge can be accomplished by using the standard Faddeev-Popov method, which introduces the 
Faddeev-Popov determinant also in the integration measure 

(*|...|y) = J DAJ(A)r(A)...Y(A). (1.6) 

2. Zero-temperature Yang-Mills theory 

One can solve the Yang-Mills Schrodinger equation in perturbation theory by expanding the 
Hamiltonian and the wave functional in powers of the coupling constant g, applying standard 
Rayleigh-Schrodinger perturbation theory From the Coulomb term (L5), which is order g 2 , 
one can extract the running coupling constant and obtains the same result as in ordinary covariant 
perturbation theory within the functional integral formulation. We are interested here, however, in 
a non-perturbative solution of the Yang-Mills Schrodinger equation and for this purpose we exploit 
the variational approach using Gaussian type wave functionals. The first variational calculations in 
Coulomb gauge were performed in Ref. [Q] and later on in Ref. [Q]. Our approach [J7|] differs from 
previous variational calculations in Coulomb gauge in the ansatz of the vacuum wave functional, 
in the treatment of the Faddeev-Popov determinant (treated fully in our approach and at least par- 
tially neglected in previous approaches) and in the renormalization, for more details see Sect. IID 
of Ref. [p. 

The variational approach developed in Tubingen uses the trial ansatz for the vacuum wave 
functional [f7|] 

y(A) = - exp 
'Det(-Z)d) 



-- / Ao.)A 



(2.1) 



which contains besides the exponential the inverse square root of the Fadeev-Popov determined. 
This has the advantage that in expectation values the Faddeev-Popov determinant in the integration 



measure, see Eq. (1.6), is cancelled. Furthermore, the static gluon propagator is with this wave 



functional given by the inverse of the variational kernel CO 

(AA) = (2coy\ (2.2) 

which shows that CO has the meaning of the gluon energy. Minimizing the vacuum expectation 
value of the Hamiltonian, (y\H\\j/) — > min, one finds []8|] for the gluon energy CO the result shown 
in Fig. [T] (dashed line). At large momenta it behaves like the photon energy m(k —> °°) ~ k and 
is infrared divergent co(k — > 0) ~ l/k. This, of course, is the manifestation of the confinement of 
gluons. Figure [I] shows the lattice results [Q] for the gluon propagator (2.2), which can be fitted by 
Gribov's formula. 

w (fc) = ^2 + — (2.3) 

with the effective mass M 880 MeV. Also shown in this figure is the result of the variational 
calculation |j|] : as can be seen the gluon propagator obtained from the variational calculation agrees 
quite well with the lattice data in the infrared and the ultraviolet, while there is some missing 
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Figure 1: Comparison of the gluon propagator 
1 /(2o) with Gaussian [ph (dashed line) and non- 



Gaussian [10 1 (full line) functional to the lattice 
data 101. 



Figure 2: Form factor of the three-gluon ver- 
tex for orthogonal momenta and comparison to 
lattice data for the 3-dimensional Landau-gauge 
vertex riFfll. 



strength in the mid-momentum regime around 1 GeV. This can be traced back to the absence of 



the gluon loop, which escapes the variational calculation with the Gaussian type of ansatz (|2. ip. In 



Ref. [10] the variational approach was extended to non-Gaussian wave functionals including up to 



quartic terms in the gauge field 

\ W [A]\ 2 = exp(-5[A]), S[A] = j ' coA 2 + i J Y (3) A 3 + i| 7 (4) A 4 . (2.4) 

To capture the gluon loop in the variational calculation one needs to include at least the three- 
gluon term y^> in the exponent of the wave functional. Then one finds the static gluon propagator 
shown in Fig. [I] (full line), which gives a substantial improvement compared to the propagator 
obtained with the Gaussian trial wave functional. Figure || shows the result for the three-gluon 
vertex. Also shown are the lattice result for the three-gluon vertex calculated in 3-dimensional 



Landau gauge Yang-Mills theory [11]. This theory corresponds to the use of the approximate 



wave functional (1.3) in 3 + 1 dimensional Yang-Mills theory in the Hamiltonian approach. Since 



the expression ( |1.3| ) represents a good approximation to the true Yang-Mills wave functional in the 



low-momentum regime (see Ref. Q12|]) we expect the 3-dimensional Landau gauge lattice result to 
agree well with the static three-gluon vertex in D = 3 + 1 Yang-Mills theory. Indeed we find a 
quite reasonable agreement between the result of the variational calculation and the lattice data in 
the infrared. 

The above represented solution correspond to the so-called critical (or scaling) solution, where 
the horizon condition d~ l (k = 0) = was assumed for the ghost form factor d(k) defined by the 
ghost propagator 

((-Ddy 1 )=d{k)/k 2 . (2.5) 

Figure || shows the result of the variational calculation for the gluon energy (shown already in 
Fig. II]) and the ghost form factor of the full variational calculation and the one with the Coulomb 



term of the Hamiltonian, Eq. (1.5), excluded. One observes that the effect of the Coulomb term is 



very small. We will therefore neglect this term in subsequent numerical calculations. 
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3. Finite-temperature Yang-Mills theory 

To extend Yang-Mills theory to finite temperatures [13, [14[] we consider the grand canonical 
ensemble with vanishing chemical potential, which is defined by the density matrix 

D = exp[-H/(k B T)], (3.1) 

where k s is Boltzmann constant. For the evaluation of the thermal averages 

<...)r = 5^ (3-2) 
we need a complete basis of the gluonic Fock space, which we choose in analogy to the ground 



state wave functional (2.1) in the form 

\k)= i 1 =\k), (3.3) 
yDet(-Dd) 

where the \k) form a complete set of states of the gluonic Fock space, which we choose in the 
following way: we decompose the gauge field (in momentum space) in terms of creation and 
annihilation operators 

A(*) = -=L=( a (*) + flt(-*)), (3.4) 

y/2(Q(k) 

where co(k) is an arbitrary (positive definite) kernel. Defining the vacuum state \k = 0) by 

a(Jfc)|0)=0 (3.5) 
this state becomes in coordinate representation 

(A|0}=exp(-iyAfflAV (3.6) 

A complete basis of the gluonic Fock space is then given by 

|0>, a f (k)\0), a\k)a\k')\0), ... (3.7) 
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The exact density matrix ( f3~l] ) is too difficult to handle given the complicated form of the Yang- 
Mills Hamiltonian. For this purpose we replace the Yang-Mills Hamiltonian in the density matrix 
by a single-particle Hamiltonian 

D = exp[-h/(k B T)], h = J &kQ.(k)a\k)a{k), dk = d 3 k/(2n) 3 . (3.8) 



Since h is a single-particle Hamiltonian the thermal averages (3.2) with the density matrix (3.8) can 



be worked out by using Wick's theorem. For the gluonic occupation numbers one finds 

(a(k)a* (k')) T = {Inf 8{k-k')n{k), with n(k) = [exp(£2(jfc)/(jfc B r)) - l] _1 , (3.9) 

which are the usual Bose occupation numbers with Cl(k) representing the single-particle energies. 
By means of Wick's theorem all thermal averages can then be expressed in terms of the gluon 
propagator 

(A(k)A(k')) T = (2n) 3 8(k + k , )(l+2n(k))/(2(o{k)). (3.10) 
From the density matrix ( |3~8| ) we find the entropy S and the free energy F 

S = — & B TrDlnD, F = {H) T -TS. (3.11) 



So far the two kernels Q.(k), entering the density matrix (3.8), and co(k), entering our basis states 



(3.6), are completely arbitrary. We now determine Cl(k) by minimizing the free energy (3.11). 
Instead of taking the variation with respect to Cl(k) it is more convenient to take the variation with 
respect to the finite temperature occupation numbers n(k) ( |3.9| ), which is equivalent since n(k) is a 
monotonic function of Q.(k). This yields 

5e[n,co] 

m= ^w> (3 - 12) 

where e[n,(o] = (H) T /(2(N 2 - 1) -V) is the energy density per gluonic degree of freedom (V is 
the spatial volume). From the analogy with Landau's liquid Fermi theory we identify Q.(k) as a 
quasi-gluon energy. Evaluating the thermal expectation value of the Yang-Mills Hamiltonian {H)j 
up to two loops one finds for the quasi-gluon energy 

The kernel co(k) can be chosen, in principle, completely arbitrary (except that it has to be 
positive definite) and the results of our variational calculation should not depend on the choice of 
(o{k). However, since we have introduced approximations (calculating the energy up to two loops) 



our results do depend on co{k). The best we can do is to vary the free energy F Eq. ( 3.1 1 ) with 
respect to co(k). This guarantees that the free energy is at least to first order independent of co(k). 
From 8F / 8co(k) = we find the gap equation 

(0 2 (k)=k 2 + x 2 (k)+I tad [n]+I c [n](k), (3.14) 

where #(&) is the ghost loop, 7 t adH is the tadpole and I c [n](k) the one gluon-loop contribution 
from the Coulomb term, see Fig. ^ The finite-temperature modifications arise exclusively from 
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(a) (b) (c) 

Figure 4: Diagrammatic representation of the contributions to the gap equation: (a) full ghost loop, (b) 
tadpole, (c) one-gluon loop contribution from the Coulomb term. Small black dots and open circles represent, 
respectively, bare and full vertices of the Hamiltonian; the double line denotes the Coulomb kernel. 




Figure 5: Diagrammatic representation of the ghost DSE. 



the finite-temperature part of the gluon propagator, Eq. ( 3.10 ), which depends on the occupation 
number n(k). The finite-temperature modifications are all ultraviolet finite so the renormalization 
of the gap and Dyson-Schwinger equations can be done in exactly the same way as at zero tem- 
perature, see Ref. [15]. The gap equation has to be solved together with the Dyson-Schwinger 
equation for the ghost propagator, which is illustrated in Fig. |5[ This equation is the same one as at 
zero temperature except that the gluon propagator is replaced by its finite-temperature counterpart 
( 3.10 ). As we have illustrated above for zero temperature, the Coulomb term of the Yang-Mills 
Hamiltonian does barely influence the ghost and gluon propagator. Therefore we will ignore the 
Coulomb term in the following. 1 Neglecting the Coulomb term leads to substantial simplifications 
[14]. The quasi-gluon energy H(jt) of the density matrix (3.8) and the kernel (o(k) of the vacuum 
wave functional (3.6) become then identical. Furthermore, the ghost Dyson-Schwinger equations 
and the gap equation then decouple from the Dyson-Schwinger equation for the Coulomb form 
factor [§]. 



4. Infrared analysis 



Before we present the numerical results let us summarize the results of an infrared analysis for 
the remaining gap and ghost Dyson-Schwinger equations. Let us first recapitulate the result of the 



infrared analysis at zero temperatures [16]. 



For the infrared analysis we assume power-law behaviours of the propagators involved 

co(p->0) =A/p a , d(p^0)=B/pP. (4.1) 

Assuming the horizon condition d~ Y (p = 0) = 0, which implies /3 > 0, we find from the ghost DSE 
the sum rule 

a = 2[5+2-d, (4.2) 



'We should, however, stress that the Coulomb term is utterly important for the quark sector since it provides the 
confining potential for the quarks, see further below. 
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where d is the number of spatial dimensions. For Coulomb gauge in d = 3 we obtain a = 2/3 — 1, 
and including also the gap equation one finds the following two solutions 

/3 = 1 (1.001), J8 = 0.796 (0.794), (4.3) 

where the results from the numerical evaluation are given in the bracket. In d = 2 dimensions one 
finds a single solution, which within the angular approximation is given by 

P = 0.5 (0.45) (4.4) 

while without the angular approximation one finds /3 = 0.4. 

At arbitrary finite temperature the infrared analysis cannot be carried out since the gluon en- 



ergy Q(k) = co(k) enters the thermal occupation numbers n(k), Eq. (3^9), in the exponent. However, 



at very high temperatures we can expand this exponent and the occupation numbers reduce to 

n{k) = k^T/(o{k). (4.5) 
Then co{k) enters only as power in the gap and DSE and we can carry out the infrared analysis in 



the standard fashion. One finds the same sum rule Eq. (4.2) as at zero temperature. However, one 
obtains only a single solution for the infrared exponent 

= 1/2, a = 0. (4.6) 

This is the infrared exponent for the ghost in d = 2 dimensions, which reflects the fact that at high 
temperature a dimensional reduction to the 2 + 1-dimensional theory occurs. 



5. Numerical results 

Figure |6] shows the ghost infrared exponent /3 determined from the numerical solutions as 
function of the temperature. As one observes the infrared exponent stays constant below a critical 
temperature where it suddenly drops and approaches the value /3 = 1/2 for high temperatures, 
in agreement with the infrared analysis. Both zero-temperature solutions, j8 = 1 and /3 = 0.796, 
convert to the same high-temperature solution. Figures [7] and |8] show the numerical solutions for 
the ghost form factor and the gluon energy for zero temperature and above the deconfinement 
temperature. One observes that the infrared behaviour is indeed in accord with the findings of the 
infrared analysis. The sudden drop of the infrared exponent of the ghost form factor (see Fig. |6|) can 
be used to define the deconfinement phase transition temperature T c . Fitting at T = the numerical 
solution for (o(k) to the lattice gluon propagator [Q] to fix the physical scale one finds for the critical 
temperature of the deconfinement phase transition 

r c ~270...290Mev, (5.1) 

for SU (2), which compares well with the lattice result of T c = 295 MeV. 

Figure ^ shows the infrared mass of the gluon defined by the gluon energy at the numerical 
infrared scale Air, 

m(T) = co(k = X 1R ,T). (5.2) 
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Figure 6: IR exponent of ghost as function of 
the temperature. 
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Figure 7: The ghost form factor. 
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Figure 8: The gluon energy. 



Figure 9: The IR mass of the gluon. 



This mass behaves similarly as the infrared exponent of the ghost form factor. It stays constant 
below a critical temperature, where it suddenly drops and after that starts rising linearly with the 
temperature. Zooming into the behaviour of m(T) in the transition regime of Fig. ^, see Fig. 1C, 
one finds for the critical exponent of the effective gluon mass defined by 



m(r)~(r/r c -i) K 



(5.3) 



a value of K ~ 0.37. This compares well with the result of K = 0.41 obtained in Ref. [17], where 
a quasi-gluon picture has been used to fit the lattice results for the energy density and the pressure, 
and assuming furthermore input from the d = 3 Ising model, which is in the same universality class 
as SU (2) gauge theory. 

Figure [ll] shows the running coupling constant at zero temperature and above the deconfine- 
ment phase transition. One observes a substantial reduction of the low energy plateau. 
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Figure 10: Critical behaviour of the effective Figure 11: Running coupling constant from the 
gluon mass. ghost-gluon vertex. 

6. Hamiltonian approch to QCD 



When the quarks are included the Hamiltonian ( L4 ) has to be supplemented by the quark term 



q\x) 



alp + gA) +f3m q(x 



(6.1) 



where q(x) denotes the quark field operator and mo is the current quark mass. Furthermore, the 
matter charge density in the Coulomb Hamiltonian (1.5) becomes p,^(x) = (x)t a q(x), where t a 



are the generators of the gauge group in the fundamental representation. For the quark sector we 



use the following ansatz for the vacuum wave functional Q18Q 



I0)c 



exp 



q f (Sp+Va-A)q 



|0>, 



(6.2) 



where |0) is the perturbative quark vacuum state, which describes a system of free quarks with 
mass mo. Furthermore S and V are variational kernels to be determined by minimizing the energy 
density. For V = Eq. (6^) defines a BCS-type wave functional considered in Ref. [|^]. The new 
important aspect of the wave function ( |6.2[ ) is the coupling of the quarks to the gauge field with 
form factor V. Without this term the quark-gluon coupling in the quark Hamiltonian ( |0| ) escapes 
the expectation value. The total QCD wave functional is then given by 



(6.3) 



IV)ym®|0)q> 

where \y)ym is the Yang-Mills vacuum wave functional given in Eq. ( |2.1| ). We keep CO fixed at 
its form determined from the Yang-Mills sector. We minimize the quark energy ((H C1 ) Q } YM with 
respect to the kernels S and V. The result of this variation is shown in Fig. 12. Figure |13| compares 
the scalar form factor with and without taking into account the quark-gluon coupling. With the 
quark-gluon coupling switched off (V = 0) we find the value for the quark condensate 

(qq) = (-113 MeVVac/aw) 3 , (6-4) 



which corresponds to the result of Ref. [19], while with the quark gluon coupling included (V ^ 0) 
we obtain 

(qq) = (-135 MeV v^c/aw) • (6-5) 
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Figure 12: The variational kernels S and V. 



Figure 13: Comparison of scalar kernel S(k) with 
(V ^ 0) and without (V = 0) the quark-gluon cou- 
pling included. 
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Figure 14: Running quark mass. 



This is an increase of the figure in the bracket by about 20%. In Eqs. ( p.4[ ) and ( p.5| ) a c and a w 
denote respectively the Coulomb and Wilsonian string tension. Lattice results show that the ratio 
of these quantities is given in the range of 



Oc/Ow ^2...3. 
This implies a quark condensate in the range of 

(qq) = (-191... 234 MeV) 3 , 
which compares well with the phenomenological value of 

{qq) = (-230 MeV) 3 . 



(6.6) 



(6.7) 



(6.8) 



Figure [14] shows the dynamical quark mass as a function of the momentum with and without 
the quark gluon coupling. We observe an essential increase in the quark mass when the coupling of 
the quarks to the gluons is included. For the constituent quark mass defined by the zero momentum 
value of the running quark mass M = M(k = 0) we obtain without the quark-gluon coupling (V = 0) 



M = 84 McVa/ocT^ 



(6.9) 
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and with the quark-gluon coupling included (V ^ 0) 



M= 132 MeVA/cTc/aw, (6.10) 



which is an increase by 57%. Assuming again the lattice result Eq. ( p.q ) for the ratio of the string 
tensions we obtain Jl8| ] constituent masses in the range of 

M= 187...230MeV, (6.11) 

which brings the constituent mass into the region of its phenomenological value. 

The results obtained so far in the variational approach to QCD in Coulomb gauge are very 
encouraging and call for more detailed investigations. 
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